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Abstract. We describe XMM-Newton Guaranteed Time observations of a sample of eight high redshift (0.45<z<0.62) clusters. 
The goal of these observations was to measure the luminosity and the temperature of the clusters to a precision of ~10%, 
leading to constraints on the possible evolution of the luminosity-temperature (Lx - T^) relation, and ultimately on the values 
of the matter density, Qu, and, to a lesser extent, the cosmological constant Qa- The clusters were drawn from the SHARC 
and 160 Square Degree (160SD) ROSAT surveys and span a bolometric (0.0-20 keV) luminosity range of 2.0 to 14.4 xlO'*'' 
erg s ' (Ho=50,fiM=l> S^a=0). Here we describe our data analysis techniques and present, for the first time with XMM- 
Newton, a Lx - T^x relation. For each of the eight clusters in the sample, we have measured total (r < ryi,-iai) bolometric 
luminosities, performed yS-model fits to the radial surface profiles and made spectral fits to a single temperature isothermal 
model. We describe data analysis techniques that pay particular attention to background mitigation. We have also estimated 
temperatures and luminosities for two known clusters (Abell 2246 and RXJ1325. 0-3814), and one new high redshift cluster 
candidate (XMMU J084701. 8+345117), that were detected off-axis. Characterizing the - T„ relation as = L^C^)", we 
find Le = 15.9!j-^ x 10'"*erg s"' and Q'=2.7±0.4 for an Qa = 0.0, Am = 1-0, H g = 50 km s"' Mpc"' cosmology at a typical 
redshift z ~ 0.55. Comparing with the low redshift study by |Markevitch, 1998] we find a to be in agreement, and assuming 
Lx - Tx to evolve as (1 -I- zY", we find A=0.68±0.26 for the same cosmology and A = 1.52;^|5 2^ for an Oa = 0.7, Qm = 0.3 
cosmology. Our A values are very similar to those found previously by |Vikhlinin et al., 2002 'using a compilation of Chandra 
observations of 0.39 < z < 1.26 clusters. We conclude that there is now evidence from both XMM-Newton and Chandra for 
an evolutionary trend in the - relation. This evolution is significantly below the level expected from the predictions of 
the self-similar model for an Ha = 0.0, Dm = 10, cosmology, but consistent with self-similar model in an ^a = 0.7, = 0.3 
cosmology. Our observations lend support to the robustness and completeness of the SHARC and 160SD surveys. 

Key words. X-rays:Galaxies: clusters : 



1. Introduction 

1.1. Motivation 

As the most massive gravitationally bound objects in the 
universe, galaxy clusters are particularly sensitive to the 
evolution of the density perturbations responsible for their for- 
mation. Cluster abundance as a function of mass and redshift 
is dictated by the mass function ( |Press and Schech ter, 1974| 
[Jenkins et al., 2001} , which gives the comoving space density 
of collapsed objects as a function of mass and redshift. In 
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standard models, the mass function falls off as a Gaussian 
at the high mass end, reflecting the Gaussian nature of the 
density perturbations. Galaxy cluster abundance is therefore 
exquistely sensitive to the amplitude of the perturbations 
and its evolution. Since this evolution is controlled by the 
underlying cosmological background, observations of clus- 
ter abundance offer an effective way to constrain certain 
cosmological parameters, such as the density parameter 
( lOukbir & Blanchard, I992|BIanchard and Bartlett,T998i l 
or the dark energy parameter and equation-of-state 
( Majumdar & Mohr, 2003 Wang & Steinhardt, 1998 ). The 
present day cluster abundance is degenerate in the matter 
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Fig. 1. The predicted X-ray temperature function at redshifts 
z = 0.05 (black), 0.33 (light grey) and 0.5 (dark grey) for a flat 
high density model (left; = 0.12) and a concordance model 
(right; Q.^ = 0.7), both fitted to the local (z ~ 0.05) temperature 
function. 



density and amplitude of the power spectrum; evolution of 
the abundance breaks this degeneracy. Constraints obtained in 
this manner are complementary to others that essentially rely 
on determinations of cosmological distances; for example, 
those found by observations of supernovae type la or by 
measurements of cosmic microwave background anisotropies. 

Cluster mass is, however, difficult to measure directly, and 
in practice one seeks a direct observable that is closely re- 
lated to virial mass. Lensing surveys would seem the most 
suited to the task, as the efi'ects of lensing are of course di- 
rectly related to mass (albeit projected along the line-of-sight). 
However, cluster mass estimates from weak and strong lens- 
ing remain controversial as they suffer from several system- 
atic uncertainties (e.g. projection effects). Among X-ray ob- 
servables, intracluster gas temperature is expected to be tightly 
correlated with virial mass, an expectation borne out by sim- 
ple hydrostatic considerations as well as by numerical simula- 
tions dEvrard et al., 1996||Bryan and Norman,"?998t . X -ray lu- 
minosity, on the other hand, is a much less robust mass indica- 
tor, despite being significantly easier to measure, because it de- 
pends on the density profile of the intracluster gas, the physics 
of which is currently difficult to model. 

With a calibrated - M (temperature-mass) relation, 
the mass function can be translated into an observable clus- 
ter temperature distribution function dn/dT. The exact T - 
M relation to use is of course a key ingredient, one that 
may be addressed, for example, using numerical simulations, 
or directly from detailed observations that determine both 
cluster mass and temperature (e.g. Nevalainen et al., 2000j . 
Figure [2 compares temperature function predictions for a 
high matter density and a flat model, both normalized to the 
present-day, obse rved dn/dT (e.g. [Hen ry and Arnaud, 1991] 
[Edge et al., 1990| [Blanchard et al., 2000| |Ikebe et al., 2002t ' 
Evolution toward higher redshift is strikingly different in two 



cosmological models differing by their matter content alone 
and illustrating the above argument. 

Despite this obvious promise, it has been difficult to use 
dn/dT(z) to explore the cosmological model because direct 
measurements of Tx at z > require long integration with satel- 
lite observations. Prior to the launch of Chandra and XMM- 
Newton, the number of measured cluster temperatures at high 
redshifts was very small. This is attested to by the fact that 
the most distant determination of dn/dT ( |Henry, 2000| to date 
was based on a sample of only fourteen z > 0.3 clusters (mean 
~0.38) from the Einstein Medium Sensitivity Survey (EMSS; 
IGioia et al., 1990/ 

An alternative, but related approach is to apply a 
luminosity-temperature (Lx - T^) relation to a flux- 
limited sample, thereby obtaining either the temperature 
function, or a redshift distribution at given temperature 
dOukbir & Blanchard, 199 7 Sadat et al.,1998 Reichart et a l., 1999|Borgani et 
The advantage of this approach is that it does not require de- 
tailed X-ray spectroscopy of all the clusters in a flux-limited 
sample to establish the ~ T^x relation and its evolution. 
There are now several flux limited catalogs of medium 
to high redshift clusters, in addition to the EMSS, to 
which one can apply this technique. These include those 
based on data from the ROSAT All Sky Survey (e.g. 
lEbeling, Edge, & Henry, 2001' 'Zhang et al.,"!^^ ), NEP (e.g. 
Gioia e t al.(2001 ) He nry et al.(2001)i; and from the ROSAT 
pointing archive (e.g. |Romer et al., 2000| [Burke et al., 2003| 
iVikhlinin et al., 1998| [MuUis et al., 20031 IScharfetal. 19971 
[Rosati et al.,"T995} . 



1.2. The XMM-Newton a-Project 



The goal of the XMM-Newton Q-project (Bartle tt et al., 200U 
is to increase the number of high quality X-ray cluster tem- 
perature measurements at z > 0.3, and thus enhance the sci- 
entific yield from the various ROSAT surveys. We describe 
below XMM-Newton observations of eight ROSAT clusters 
at z > 0.4 (median of z = 0.54) performed as part of the 
Guaranteed Time programme. Related observations conducted 
during the open time programmes of seven 0.3 < z < 0.4 
clusters are discussed elsewhere ( [Majerowicz, et al., 2002(b)'l 
'Majerowicz, et al., 2004). Cosmological interpretation is pre- 
sented in Vauclair et al, 2003. In Sections 2 & 3 we describe the 
observations and data analysis techniques. In Section 4 we dis- 
cuss each of the eight clusters in turn. In Section 5 we present 
our Lx - T^x relation and compare it to previous work. 

This paper serves also to introduce the Project, and repre- 
sents an opportunity to provide detailed descriptions of analysis 
techniques used in the Project and which can be be used gener- 
ally for XMM-Newton observations of cluster targets. Except 
where explicitly stated elsewhere we use a Q.m = 1, Oa = 0, 
qo = 0.5 model with Ho - 50 km s"' Mpc"' that has been 
most frequently used in the past as the parameter set for X-ray 
cluster studies. When necessary for examining cosmological 
implications we correct our results to a concordance model. 
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2. Observation Programme 

2.1. The Sample 

Table summarizes the locations, date of observations and 
other details of the eight clusters in this program. Seven of 
the eight clusters were drawn from the SHARC Surveys, 
four from the Southern SHARC (Burke et al., 20 03 1 sample 
and three from the Bright SHARC (Ro meretal., 20 00 1 sam- 
ple. The SHARC cluster samples are based on searches for 
clusters serendipitously detected in ROSAT PSPC observa- 
tions. They complement each other in that they cover, re- 
spectively, 17.7 degrees^ to a flux limit of - 3.9 x 10"'^ 
erg s 'cm^^ ( |Comns et al., 19971 and 178.6 degrees^ to a 
flux limit of ^ 1.4 X 10 erg s 'cm^^ jNichol et al., 1 999 1. 
This strategy has yielded a combined cluster catalogue that 
straddles L* over the redshift range 0.2 < z < 0.8 and 
shows a consistent picture of a non-evolving cluster lumi- 
nosity function dCoUins et al., 1997| [Burke at al., 1997 | ex- 
cept, possibly, at luminosities greater than = 5 x 
lO'*'^ erg s-i (Nic hol et al., 1 999 Ad arni et al., 2000l. Two as- 
pects of the SHARC surveys, makes them particularly well 
suited to the XMM-Newton Q project; they have been sub- 
jected both to detailed optical follow-up CRo mer et al., 2000| 
IBurke et al., 2003t and to extensive selection function sim- 
ulations jAdami et al., 2000| [Burke et al., 2003t . The eighth 
cluster in our observation programme was taken from 
flie 160 Square Degree ROSAT Survey (160SD hereafter, 
IVikhUnin et al., 1998| and |Mullis et al., 2003t . There is consid- 
erable overlap, both in terms of methodology and cluster mem- 
bers, between the 160SD and the SHARC surveys (3 of the 7 
SHARC clusters are also members of the 160SD catalogue), 
however, this particular cluster (RXJ0847.2) was not a member 
of either SHARC sample because its host PSPC observation 
did not meet the SHARC exposure time criterion. The selec- 
tion of clusters sampled in this programme was driven solely 
by a requirement to observe all z>0.5 targets from the SHARC 
surveys. Following the visibility and observability constraints 
of the XMM-Newton Guaranteed Time programme we added 
two clusters with slightly lower redshift from the SHARC sam- 
ple in addition to the one from the 160SD catalogue. 

2.2. XMM-Newton 

XMM-Newton ( |Jansen et al., 2001) comprises 3 co-aligned 
telescopes, each with effective area at l.SkeV of ~1500cm^, 
and Full Width Half Maximum (FWHM) angular resolution 
of ~5 arc-seconds. This combination of the highest ever fo- 
cused X-ray collection area, and the ability to resolve clusters at 
all redshifts, makes XMM-Newton the best suited observatory 
for this programme. The 3 telescopes each have a focal plane 
CCD imaging spectrometer camera provided by the EPIC con- 
sortium. Two also have a reflection grating array, which splits 
off half the light, to provide simultaneous high resolution dis- 
persive spectra. These two telescopes are equipped with EPIC 
MOS cameras ( |Tumer et al., 2001 J , which are conventional 
CMOS CCD-based imagers enhanced for X-ray sensitivity. 
The third employs the EPIC PN camera ( IStriider etal., 200U 



which is based on a pn-junction multi-linear readout CCD. The 
EPIC cameras offer a field of view (FOV) of ~30 arc-minute di- 
ameter, and an energy resolution of typically 100 eV (FWHM) 
in the range -0.2-10 keV. The 7 CCDs in each of the MOS 
cameras are about 10 arc-minutes square each. The central chip 
encompasses the whole of ourz > 0.4 clusters. The 12 CCDs in 
the PN are about 4x13 arc-minutes. Even if correctly centered 
on the boresight PN CCD, some portion of the target cluster 
emission may spill onto neighbouring CCDs, and across dead 
zones between CCDs. In each camera, an aluminized optical 
blocking filter was deployed, the thickness of which was cho- 
sen to suit the expected brightness of nearby serendipitous ob- 
jects in each field. Seven of the eight clusters were observed 
in the Full Frame Imaging mode, as appropriate for weak ex- 
tended targets, but one (RXJ1325.5) was observed when the 
MOS cameras were in "Window" mode, due to the presence of 
an unusually bright nearby point source. 

The data were processed through the XMM-Newton 
Science Analysis Sub-System (SAS; [Watson et al., 2001t ver- 
sion 5.3, in order to register photons from detector to sky co- 
ordinates, to correct energy data for gain and charge transfer 
losses, and remove instrument noise artifacts. This provided 
calibrated event lists as a starting point for our detailed data 
reduction. 

3. Data Reduction 

3.1. Rate Filtering 

Part of the XMM-Newton orbit lies within the magnetosphere, 
and consequently the spacecraft can encounter clouds of pro- 
tons accelerated by magnetic reconnection. When these par- 
ticles scatter through the mirror system they are concentrated 
onto the focal plane, and an enhanced background rate can oc- 
cur These intervals were identified by forming histograms of 
events with energy >10 keV, located in single pixels, in time 
bins of 50 (100) seconds in the PN (MOS) camera(s). Due 
to variations in the baseline raw cosmic ray rate experienced 
through the mission to date (probably due to modulation by the 
solar activity and/or the seasonal variation of satellite apogee 
direction), we prefer not filter at fixed background count rates, 
rather for each exposure we define +3cr limits after ignoring 
the highest count rate periods (Figure |2 and |3- Table |5] sum- 
marizes the total on-axis Good Time Interval (GTI) exposure 
times extracted for each observation using this procedure. 

For clusters RXJ 1354.2-0222 and RXJ1701. 3+6414 the 
high prevalence of proton flares mitigated against the auto- 
mated filtering procedure, and a manual analysis of the the 
"curve of growth" for signal:noise in the cluster target region 
was employed. 

3.2. Vignetting 

The clusters in our sample are expected to be extended on 
a scale of 2-3 arc-minutes, and in order to assess accurately 
the surface brightness distribution, a small correction for en- 
ergy dependent telescope vignetting must be applied. This was 
performed by the SAS task EVIGWEIGHT which assigns a 
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Cluster ID 


RA 


Dec 


2 


Date 


Obs ID 


Duration 


Piltcr 




n3-37-4S 


Z.^ . Z.VJ 


0.577' 


7001-08-1 8T1 1 •46-53 


01078(S0401 


58942 


MFDTT TM 


RXJ0505. 3-2849 


05:05:20 


-28:49:05 


0.509' 


2001-09- 01X13:17:36 


0111160201 


48867 


XHIN 


RXJ0847.2+3449 


08:47:11 


+34:49:16 


0.560'' 


2001-10- 07X11:55:10 


0107860501 


81143 


XHIN 


RXJ 11 20. 1+43 18 


11:20:07 


+43:18:05 


0.600'' 


2001-05-08X20:50:37 


0107860201 


22627 


XHIN 


RXJ1325.5-3826 


13:25:20 


-38:24:55 


0.445' 


2002-01-19X02:30:04 


0110890101 


60894 


MEDIUM 


RXJ1334.3+5030 


13:34:20 


+50:30:54 


0.620* 


2001-06-07X20:19:43 


0111160101 


47614 


XHIN 


RXJ1354.2-0222 


13:54:17 


-02:21:46 


0.551' 


2002-07- 19X15:16:53 


0112250101 


33374 


XHIN 


RXJ1701. 3+6414 


17:01:23 


+64:14:08 


0.453* 


2002-05- 31X17:49:42 


0107860301 


18172 


MEDIUM 



Table 1. Summary of pointing directions for each observation, the overall scheduled durations, and redshift data (" Southern 
SHARC, .Burke et al., 2003. * Bright SHARC, .Romer et al., 2000. " 160SD ,VikhUnin et al., 1998) . 
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Fig. 2. Count rate in PN (top) and MOS (bottom) cameras 
in the RXJ 1120.1 observation after selection for single pixel 
events above lOkeV. Typical rates ~15cts/50 sec bin (PN) and 
~8cts/100 sec bin (MOS) 
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Fig. 3. Histogram of the count rate sequence for PN (top) and 
MOS (bottom) cameras in the RXJ 11 20.1 observation after se- 
lection for single pixel events above lOkeV. The dashed lines 
indicate the upper and lower (+3cr) bounds used for the Good 
Time Interval filtering. 



Cluster ID 


MOS 


PN 


% lost 


RXJ0337.7-2522 


23421 


23611 


60 


RXJ0505. 3-2849 


31851 


24983 


35 


RXJ0847.2+3449 


43793 




50 


RXJ1120.1+4318 


18333 


16208 


18 


RXJ1325. 5-3826 




48753 


20 


RXJ1334.3+5030 


40713 


34909 


15 


RXJ1354.2-0222 


9090 


7248 


73 


RXJ1701. 3+6414 


9379 


4382 


50 



Table 2. Summary of the usable exposure durations obtained 
for the different clusters after accounting for the GTI filtering, 
telemetry losses etc.. In some observations data from only one 
camera was available. For operational reasons the PN camera 
commences exposures somewhat later than the MOS cameras. 
Losses due to proton flares, telemetry drops etc. are given in 
column 4 as a percentage of the total exposure time. 



weighting value to each detected photon, accounting for tele- 
scope and CCD efficiency variations (the latter being negligi- 
ble by comparison). Thereafter, spectrum extraction and image 
product generation is automatically weighted for the vignetting 
correction, and response matrices for the on-axis location can 
be used. ' 

An exposure map was also generated for each camera ob- 
servation. These maps accounted for spatial variations, such as 
dead pixels, noisy readout columns and chip gaps, not included 
in the vignetting correction. The XMM-Newton satellite is usu- 
ally very stable during observations, so the inter-CCD gaps re- 
main approximately fixed in sky projection. It is therefore im- 
portant to coiTect for flux lost from extended cluster regions 
using the exposure maps. This is especially true for the PN im- 
ages, because the PN CCDs are smaller and the gaps closer to 
the target than those in the MOS cameras. 

3.3. Background subtraction techniques 

Background subtraction is a very important step in our analy- 
sis. Even after the proton flare removal described above (sec- 



' In principle this may not apply correctly for objects /ar from the 
CCD array centre, as the detector response redistribution matrices may 
diverge slightly from the on-axis case. It is probably also true that 
mirror vignetting calibrations become less secure with increasing off- 
axis angle. 
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tion B.U . and after masking out point sources, our cluster ob- 
servations will be contaminated, to varying degrees, by three 
different background signals; (/) the cosmic X-ray background, 
(;;) a particle background induced by incident cosmic rays and 
(///) some residual soft proton contamination. We will refer to 
these as the cosmic, particle and proton backgrounds respec- 
tively, hereafter There are two approaches to background sub- 
traction, either one can use the source observation itself (this is 
known as in-field subtraction), or one can use background tem- 
plate files. The templates are generated by combining several 
deep, blank field, observations and are especially useful when 
the source of interest covers a large fraction of the field of view 
(e.g. a nearby cluster). We have used a combination of both 
techniques in the analysis below. 

For each of the eight clusters in our study, we applied the 
same rate filtering, vignetting correction and detector-to-sky 
conversion used for the reduction of the cluster observation, 
to the background template files provided by the XMM-SOC 
( |Lumb et al., 2002 1. This allowed us to extract identical phys- 
ical detector regions for both the cluster and the background 
using the SAS task ATTCALC. We then re-normalized the 
background template to account for any differences in the par- 
ticle background count rates between the cluster and back- 
ground template observations. This step was necessary be- 
cause the particle background varies with observation epoch 
(Section f3.3.U . To make this correction, we took advantage of 
the fact that high energy particles are not focused by the tele- 
scope optics, so the particle background count rate can be mea- 
sured from areas outside the telescope FOV. 

We used re-normalized background templates (energy 
range 0.3^.5 keV) during the spatial analysis (Section l3.5> . 
We also used background templates (energy range 0.3-10 keV) 
during the spectral analysis of two of our two brightest clus- 
ters (RXJl 120.1 and RXJ1334.3). For this, we employed the 
so-called "double subtraction" technique ( Arnaud et al., 20011 . 
which involves making an additional correction to compensate 
for the fact that the Galactic Halo and the Local Hot Bubble 
component of the diffuse cosmic X-ray background varies sig- 
nificantly across the sky. To determine this correction, we com- 
pared off-axis, source free, regions in both the cluster fields 
and their corresponding background templates. We found the 
correction to be small, due to the lack of soft X-ray emission 
features at the high Galactic latitude locations of RXJl 120.1 
and RXJ1334.3. We have used the results from the double sub- 
traction analysis of RXJl 120.1 and RXJ1334.3 to validate the 
results from the in-field background subtraction, see below. 

For a variety of reasons, e.g. the choice of filter^, low signal 
to noise, non-zero off-axis angle etc., we were not able to apply 
the double subtraction technique to all our clusters. Instead we 
have to rely on the in-field technique. The radial dependence of 
the vignetting means that this technique is most successful for 
point sources, however it should still work well for our clusters 
since, at z > 0.4, they are barely more extended than the instru- 
ment point spread function (PSF). We illustrate this by consid- 

^ The |Lumb et al., 2002| background templates are only appropriate 
for observations made using the Thin filter; three of our clusters were 
observed through the Medium filter (TableQ 



ering a cluster described by a j8-model with yS -0.67 and a core 
radius of 25" (~ 185/;g,5 kpc typical values for the clusters in 
our sample, Table|5j- At a radius enclosing 75% of the counts 
from such a cluster, the vignetting is different by less than 1% 
from the aim-point value. Therefore we posit that an accurate 
representation of the cosmic X-ray background spectrum at the 
cluster position can be derived using a nearby in-field back- 
ground aperture. For the spectral analysis of our clusters we 
generally used an annular background aperture centered on the 
source. The annuli were chosen on a case by case basis, but 
were typically 2-3' wide with an inner radius no less than 
3' from the cluster center (to ensure that no flux from outlying 
regions of the cluster was erroneously removed). We note that 
in the case of RXJ1325.5, we could not use an annulus, due 
to the proximity of a very bright point source (Section r4.6> . In 
Sections [4. ll and R"^ we compare the results of the spectral fit- 
ting for the clusters where both in field and double subtraction 
could be applied. As we did not find significant differences, 
we conclude that the in field background subtraction technique 
should be valid for our sample. 

3.3.1. Impact of Residual Proton Contamination 

The ratio of the proton to particle backgrounds is not con- 
stant; even after proton flare mitigation (section lTll . the proton 
background may have a significantly different counting rate in 
the cluster observation compared to background template files 
(Mar kevitch, 2QQ7\ . This is illustrated in Figure |3 The trian- 
gles indicate the ratio of the particle background in the cluster 
observation to that in the background template file (column 2 in 
Table 13 • The squares show the equivalent information for the 
proton background (column 3 in Table|3- For this comparison, 
the particle background was estimated from the >10 keV count 
rate in CCD areas outside the telescope FOV, whereas the pro- 
ton background was estimated by subtracting that value from 
the > 10 keV count rate within the FOV. From the Figure|4] it is 
clear that the proton background varies with epoch and is anti- 
correlated with the particle background. This anti-correlation 
can be explained if 

- enhanced solar activity deposits more protons into the mag- 
netosphere, expanding the latter and thus shielding the cos- 
mic ray flux more efficiently. 

- seasonal variations in orbit take the spacecraft in and out 
of the magnetosphere at apogee, so that cosmic ray flux 
shielding varies inversely with the exposure to magneto- 
spheric protons. 

Although this anti-correlation tends to maintain a more con- 
stant total background rate, the result of this anti-correlation 
is that we cannot be certain the background template re- 
normalizations described above, compensate correctly for the 
temporal variations in the proton background. This suggests 
that the in-field background subtraction method may be prefer- 
able to double subtraction technique for spectral analysis. 
However, the in-field technique is complicated by the fact 
that the scattering of low energy protons at the miiTors oc- 
curs over angles somewhat larger than for the X-ray re- 
flection ( |Rasmussen et al., 1999| l. Therefore, when using task 
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Fig. 4. Comparison of the relative scaling of particle (triangle) 
and proton (square) background rates {after GTI filtering) in 
the eight cluster observations, with respect to the components 
of the background template files. Data averaged for all cameras 
in observation. 



EVIGWEIGHT (Section [T^, there is the potential to over- 
weight the proton background. The proton background spec- 
trum is hard, so this over-weighting could result in an artificial 
softening of the cluster spectrum and hence a lowering the es- 
timated cluster temperature, which might in turn mimic evolu- 
tion in the measured - relation (see section lsTl . 



To examine the likely impact of systematic errors in the 
proton vignetting correction, we used the EVIGWEIGHT vi- 
gnetting corrected background template files to create an image 
that was essentially free from proton contamination. We did 
so by applying very strict count-rate filtering. We then made a 
comparison image by applying less conservative count-rate fil- 
tering to the background template files. For this we used the fil- 
tering criteria derived from the RXJ1701.3 observation (which 
suffered from unusually high proton contamination); the com- 
parison is presented in Figure|5] The two count rate images we 
derived, should be identical within the noise, except with re- 
gard to the proton background. By dividing one by the other, 
and then fitting a radial surface brightness profile, we should 
therefore be able to get an impression of how the proton back- 
ground is vignetted. The results of this test (in the 0.5-7.0 keV 
band to emphasize an effect of hard proton spectrum) are shown 
in Figure [S] We also generated the equivalent figure using the 
RXJ1354.2 filtering criteria with almost identical results. From 
this investigation, we conclude that the proton background is 
over-weighted by only a few percent by EVIGWEIGHT even 
at the edges of the field of view in observations with signif- 
icant proton contamination. We therefore chose to ignore the 
over-weighting of the proton background in our analysis, ex- 
cept in the cases of RXJ1701.3 and RXJ1354.2. For these, we 
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Ratio of the surface brightness distributions in blank 
showing proton contamination in RXJ 170 1.3 +64 14 



Cluster ID 


Outside FOV 


Inside FOV 


RXJ0337.7 -2522 


0.72 


1.82 


RXJ0505.3 - 2849 


0.70 


1.27 


RXJ0847.2 + 3449 


0.68 


2.14 


RXJ1120.1 +4318 


0.99 


1.05 


RXJ1325.5 - 3826 


1.05 




RXJ1334.3 + 5030 


1.00 


1.23 


RXJ1354.2 - 0222 


0.52 


3.44 


RXJ1701.3 + 6414 


0.48 


5.13 



Table 3. The scaling ratio of the > 10 keV count rates out- 
side (column 2) and inside (column 3) the telescope FOV in 
the background template files compared to the cluster files. The 
former ratio is an indication of the particle (cosmic ray) back- 
ground rate, while the latter indicates the level of the soft proton 
background rate. These ratios are plotted as a function of Julian 
date in Figure 0] 



applied a small (~2%) scaling to the background spectrum de- 
rived from an in-field annulus around the cluster. ^ 



3.4. Image Products 

Images for each camera were compiled in the 0.3^.5 keV 
energy band, this maximizes source count rates for cluster 
spectra that are characterized by a ~4keV temperature, and 
furthermore provides a guard against corrections at the soft- 
est and hardest energy bands that are subject to the largest 
potential Galactic and cosmic ray subtraction inaccuracy. A 
spatial binning of 4.3 arc-seconds per pixel, was employed, 
which slightly oversamples the mirror FWHM. The 3 separate 
count-rate images from the EPIC cameras were exposure cor- 
rected and co-added. In each image field, we use the SAS task 
EBOXDETECT to identify point sources via a sliding box de- 



^ We caution the reader not to use Figure|5|as a general template for 
XMM-Newton analysis, as both the spectrum and the absolute count 
rate of the proton background are likely time dependent. 
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Fig. 6. (a) Central ~ 18 x 18 arc-minutes of the RXJl 120. 1 field - smoothed by a 4"cr Gaussian. The large circle represents the 
spectral extraction radius, and the smaller circles locate the point sources which were excised. Similar images for RXJ1334.3 
(b), RXJ0337.7 (c) and RXJ0505.3 (d). The RXJ1334.3 image 15 x 15 arc-minutes; the RXJ0337.7 and RXJ0505.3 images are 
14 X 14 arc-minutes. 



tection algorithm (see Table lA.il . In Figures |6landQ we show 
the inner portion of the eight cluster images after a Gaussian 
smoothing of ~ 4" (using the FTOOL f gauss). It is clear from 
these images that the clusters are well resolved and encompass 
a variety of morphologies. Point sources masked during the 
spatial and spectral analysis (section |^ are indicated by small 
circles. The larger circles denote the apertures used to generate 
the cluster spectrum. 

3.5. Radial surface brightness profiles 

Before generating the profile, we defined the centroid of the 
cluster brightness distribution using a 2-d Gaussian fit around 
the core of the raw cluster image (0.3-4.5 keV). These cen- 
troids are given in Table |5] row 1. Next, a background correc- 
tion was applied by subtracting, pixel by pixel, the correspond- 
ing re-normalized background template (Section [3.3> . The ra- 
dial bins were chosen so that the background-subtracted counts 
per bin, in the co-added profile, were at > 3cr significance. 



A yS-model convolved with the XMM-Newton telescope PSF 
appropriate at the position of the cluster centroid was then fit 
to these profiles, using simple ;^'--minimization (background 
was fixed). For the PSF convolution, we used calibration file 
XRTn_XPSF_0004.CCF in SAS medium accuracy mode, avail- 
able from the XMM-SOC calibration ftp site. The PSF was 
constructed by co-adding the contribution from PSFs at dif- 
ferent energies using a weighting scheme appropriate for a 
4keV thermal spectrum. Fortuitously the XMM-Newton tele- 
scope PSF is a rather weak function of X-ray energy, so that 
any deviations from this default spectrum would have a neg- 
ligible impact on our fits. We were able to use the on-axis 
PSF for the convolution, except in the case of RXJl 325 .5. 
For this cluster, which was observed off-axis, the appropriate 
off-axis PSF model was generated. We note that the PSF cor- 
rection was applied separately for each camera. Fitted values 
for p and the core radius are given in Table |5] rows 8 & 
9. To convert the core radii from angles to distances for this 
Table, we assumed a spatially flat cosmology with Qm = 1 
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and Hq - 50 km s ' Mpc ' . The eight radial profiles and their 
respective best fit PSF convolved /? models are presented in 
Figures[T9landl20l 

3.6. Spectral extraction 

Cluster spectra were extracted from circular regions centred on 
the cluster centroids, with radii ranging from 90 to 145 arc- 
seconds (typically 120", see Table|5] row 13). From these re- 
gions we extracted data from the event list using the recom- 
mended pixel selections. For the PN camera we used single 
and double pixel events (event patterns 0^) with all the rec- 
ommended selection flags applied. For the MOS cameras, we 
used the 0-12 event patterns. We account for background con- 
tamination in the spectra using either the double and/or in-field 
subtraction methods (Section l?3l . 

Spectral fitting was implemented with the XSPEC 
vll package l|Arnau d, 1996J , using the MEKAL mod- 
els dMewe et'aL^986| for a thermal spectrum, modified 
with interstellar absorption fMcC ammon & Sanders, 19901 
appropriate for the Galactic column density 
( |Dickey & Lockman, 1990| l. To facilitate fitting via Chi- 
square minimization, the spectral files were re-binned to 
ensure at least 25 counts per bin in order to approximate 
Gaussian statistics. The energy range was 0.3-lOkeV. In 
general, the following on-axis response files were used dur- 
ing the fits, epn_ff20^^dY9Jilter.rmf for the PN camera and 
ml(2)Jilterv9q20t5r6MllJ5.rsp for the MOSl(2) cameras. 
Exceptions to this, e.g. during the analysis of RXJ1325.5, are 
noted below. Fitted spectral parameters are quoted with Icr 
confidence limits on one interesting parameter. The spectral fits 
were typically performed simultaneously on the PN, MOS 1(2) 
spectra. Again exceptions to this, e.g. during the analysis of 
RXJ0847.2, are noted in the relevant sections. 

Following Markevitch (1998), we have also investigated the 
impact of cooling cores on our fitted temperatures, and hence 
the measured Lx - relation, by performing spectral fits after 
excising the central region. For all clusters except RXJ1701 
and RXJ1325, a 50 h^^ kpc region was excised. For RXJ1701 
& 1325 a 120kpc region was excised (see section 4.7). The 
results from these fits are given at the beginning of Table|5] See 
Section l5.2.1l for a discussion. 

3.7. Deriving Bolometric Luminosities 

For the luminosity calculations, we adopted a physically mean- 
ingful circular aperture, with a virial radius, r,,. We derived r,, 
according the T - r,. relation of Evrard et al. (1996) for each 
cluster using the best fit temperatures from the spectral fits. We 
note that uncertainties both in the measured temperature and in 
the T - r,, relation will introduce a systematic error in the de- 
rived Lx values. However, this should be insignificant, as very 
little cluster flux falls at radii close to the virial radius. 

We integrated the background-subtracted counts inside r,. 
and then corrected for any cluster flux lost in areas under 
masked point sources or inter-chip gaps. The correction factors 
are given in row 12 of Table |5] These were derived in a model 



independent fashion; the count-rate was accumulated in 1 pixel 
wide annuli centered on the cluster. If any pixels within a spe- 
cific annulus fell in chip gaps or in the point source masked, the 
annulus count-rate was scaled up accordingly. Next, we used 
XSPEC, together with the appropriate EPIC response func- 
tions, to determine the absorbed 0.3-4.5 keV flux within r^,, 
using the best fit absorbed MEKAL model from our spectral 
analysis so that it yielded the measured r < r,, count rate (the 
MEKAL normalization factors are given in row 10 of Table|5}. 
For comparison with equivalent ROSAT data such as SHARC 
or 160SD these fluxes (converted to the normal ROSAT 0.5- 
2 keV band) are listed in row 6 of Table |5] We then set the 
hydrogen column density to zero and re-determined the flux 
inside a pseudo bolometric band of 0.0-20 ke V. From the unab- 
sorbed bolometric flux, it is trivial to determine the bolometric 
luminosity for any given cosmological model. To aid compar- 
isons with previous work, we give in, row 3 of Table |5] the Lx 
values appropriate for a D.m = 1, = 0, qo - 0.5 model with 
Ho = 50 km s ' Mpc"' that has been most frequently used in 
the past as the parameter set for X-ray cluster studies. 

Also following Markevitch (1998), we estimated the cool- 
ing flow corrected bolometric luminosities by summing counts 
in a (r,„,„ < r < r^) ring. For clusters other than RXJ1701 
and RXJ1325, r„„„ = 50 h'^ kpc, while for RXJ1701&1325 
r,„,„=120 hgy kpc (see section 4.7). To account for the non- 
cooling flow flux falling at r< r„„„, we applied a cluster spe- 
cific renormalisation parameter calculated using the best fit 
model. We note that our correction approach differs from that 
of Markevitch (1998), therein a single multiplicative factor of 
1.06 was used. 

4. Individual Clusters 

We now discuss our analysis for each cluster in turn. We have 
ordered our discussion according to the Julian date on which 
they were observed (Table[l). 

4.1. RXJ1120.1+4318 

Figure |6j a) shows the vignetting corrected, background sub- 
tracted and co-added (PN+MOS1+MOS2) 0.3-4.5 keV image 
of RXJ 1120. 1+43 18. Point sources removed during analysis 
are circled and their positions listed in Table lA.ll Figure ^] 
shows the corresponding radial surface brightness distribution 
and best fit /? model; /3 = 0.77 + 0.03, Be = 27.4 + 1.2 
arcseconds (r^ - 209 /i^g kpc, Table|5}. 

From a spectral extraction region with a radius of 145" and 
using the in-field background subtraction technique, we mea- 
sured the following temperature, metal abundance and redshift 
by fitting to the PN, MOSl and MOS2 data simultaneously; 
Lx = 5.45+[;-2^ keV, Z = 0.47 + 0.09, z = 0.60 + 0.08. We note 
the consistency between the redshift measured from the X-ray 
spectrum with the optically determined value (z - 0.60; Romer 
et al. 2000). The overaU Xv of the spectral fit was 390 / 360 de- 
grees of freedom. During the fit, the hydrogen column density 
was fixed at the Galactic value (Nfj =2.1 x 10^" atoms cm"^), 
but we note that, when left as a free parameter, its best fit value 
was very similar (Nh - 2.2 + 0.4 x 10^" atoms cm"^). 
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Fig. 7. (a) Central MOS CCD image (10 arc-min) of the RXJ0847.2 field - smoothed by a 4" Gaussian. The large circle represents 
the spectral extraction radius, and the smaller circles locate the point sources which were excised. Similar images for RXJ 1325.5 
(b), RXJ1701.3 (c) and RXJ1354.2 (d). A newly discovered candidate distant cluster (XMMU J084701. 8+3451 17) is labeled to 
the upper right of RXJ0847.2. The RXJ 1325.5 and RXJ 1325.0 Southern SHARC clusters fall in an observation centered on IRAS 
13224-3809. The long steak connected to IRAS 13224-3809 comes from the "out-of-time" events collected during readout of 
the PN camera. In the RXJ1701.3 field, the nearby Abell 2246 and QSO objects are labeled. The RXJ1325.5 image is ~ 18 x 18 
arcminutes; the RXJ1701.3 and RXJ1354.2 images are ~ 14 x 14 arcminutes. 



Of the eight clusters in our sample, this object has the great- 
est signal-to-noise ratio. We have taken advantage of this to 
perform additional spectral analyses. First we have been able 
to use the double subtraction technique to investigate how the 
adopted background subtraction method impacts the spectral 
fits. Doing so, we derive = 5.6!jSf keV, Z = 0.43 + 0.06, 
z = 0.605 + 0.08 (Figure |8}. The overall Xv of the spectral 
fit was 330 / 307 degrees of freedom. During the fit, the hy- 
drogen column density was again fixed at the Galactic value 
(Nh - 2.1 X 10^" atoms cm"^), but when left as a free param- 
eter, its best fit value was (Nh - 2.2'^'^-^ x 10^'^ atoms cm"^). 
It is clear, therefore, that the spectral fits are not significantly 
changed by the choice of background subtraction technique 
(see also section l4~2l . 



We have also investigated how the spectral fits diff'er when 
we treat the PN and MOS data separately. We find them to be in 
good agreement; Tx^pn = 5.30+ 0.6 keV and Tx.mos - 5.7± 0.8 
keV. This is encouraging since, in two other cases (RXJ1325.5 
and RXJ0847.2), we do not have access to data from all three 
cameras. We have also been able to subdivide the spectral ex- 
traction region into three radial bins, and determined a crude 
temperature profile (Figure The profile is essentially flat, 
indicating that there is not a "cooling flow" region at the clus- 
ter core (this conclusion is supported by the absence of a cen- 
tral spike in the surface brightness profile). We note that for 
the PN/MOS comparison and for the radial profile, we used the 
double subtraction technique to account for the background. 
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(lO^-^ergs-') 



5.3± 0.5keV 13.9 ± 0.8 O-^^^ogt 26.4;^^ 
5.6+;j f keV 14.4 ± 0.2 0.77 ± 0.03 27.4+ 1.2 

Table 4. Comparison of parameter fits to the RXJ 11 20. 1 obser- 
vation by Arnaud et al. (2002) (first line) and this work (second 
line). Here Lx refers to the bolometric luminosity. 

4.1 .1 . Comparison with Arnaud et al. 2002 

The RXJ1120.1 observation described here have been previ- 
ously analyzed and interpreted in an earlier paper in this se- 
ries jArnaud et al., 200"2|> . The analysis procedures developed 
in that paper have formed the basis of the analysis of the eight 
clusters described here. For example, the vignetting correction 
technique used in Arnaud et al. (2002) has been subsequently 
implemented in SAS as the EVIGWEIGHT routine mentioned 
above. However, there are certain differences in our respective 
data analysis techniques, even with regard to the RXJ 1120.1 
observation. For example, we use updated calibration informa- 
tion and an updated SAS processing version. We also used a 
different method (3-sigma clipping) for the GTI filtering. We 
have also adopted the in-field background subtraction method 
as the standard for our cluster analysis (in Arnaud et al. (2002), 
only the double subtraction technique was used). Finally, we 
note our use of the both single and double PN events, com- 
pared to the selection of single events by Arnaud et al. (2002). 
In general the revised analysis of RXJ 1 120.1 has yielded very 
similar results to Arnaud et al. (2002) and we concur that, given 
its isothermal temperature profile and the absence of signifi- 
cant substructure, RXJ 1 120.1 is consistent with being a relaxed 
cluster. It is noteworthy, however, that in Arnaud et al. (2002) 
the mean temperature values determined using only MOS data 
differed from those determined using only PN data by more 
than 1 keV; T^^mos = 5.8!'-^ keV and T^,pn = 4.5!j;-^ keV. By 
contrast, we measured temperatures that differed by only 8% 
(see above). This improvement is most likely attributable to the 
improved caUbration that has become available since Arnaud 
et al. (2002) was published. For completeness we compare in 
Table |4]the values derived for various fitted parameters in the 
two analyses. For consistency with Arnaud et al. (2002), we 
quote the mean temperature derived using the double subtrac- 
tion technique (in Table|5]we quote the value from based on the 
in field subtraction technique). 

4.2. RXJ1 334.3+5030 

Figure lUb) shows the image of the RXJ1334.3-H5030 field. 
Point sources removed during analysis are circled and their po- 
sitions hsted in Table IaTTI Figure [T9l shows the correspond- 
ing radial surface brightness distribution and best fit /3 model; 
/3 = 0.66 + 0.02, 0, = 20 + 1 arcseconds (r, = 154 ± lO/ij^ 
kpc. Table |5}- 

From a spectral extraction region with a radius of 120" and 
using the in-field background subtraction technique, we mea- 
sured the following temperature and metal abundance; = 
5.20+_lH keV, Z = 0.15 + 0.08 (Figure^. The overall Xv of 
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Fig. 8. Spectral fit and residuals for RXJ 11 20.1; black - PN, 
dark and light grey - MOS1&2 (double background subtrac- 
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Fig. 9. The spatially resolved temperature profile of RXJl 120 



the spectral fit was 464 / 473 degrees of freedom. The hydro- 
gen column was fixed at the Galactic value (N// = 1.05 x 10^" 
atoms cm"^), but we note that, when left as a free parameter, 
its best fit value was very similar (Nh = 0.8 + 0.5 x 10^" atoms 
cm"^). Likewise, the redshift was fixed at its optically derived 
value of z = 0.62 (Romer et al. 2000), but when left as a free 
parameter, its best fit value was 0.63 + 0.02. The derived values 
for the bolometric luminosity and absorbed flux inside the viral 
radius are given in Table|5l We were also able to apply the dou- 
ble subtraction technique to these data and derived a consistent 
mean temperature value; = 5.05 + 0.3 keV. 

4.3. RXJ0337.7-2522 

Figure |6jc) shows the vignetting corrected, background 
subtracted image of the RXJ0337. 7-2522 field. Two ex- 
cised sources within the spectral extraction radius, XMMU 
J033742.9-252208 & XMMU J033745.9-252206, are ten- 
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Fig. 11. Spectral fit and residuals for RXJ0337.7; black - PN, 
dark&light grey - MOS1&2 



tatively associated with two stellar objects on the USNO 
Catalogue. The former is U0600-01432100, also identified as 
blue stellar object PHL4470, the latter is identified as U0600- 
01432383 (a 16.5 r magnitude object). Figure ^] shows the 
radial surface brightness distribution and best fit /? model; 
/?=0.76+[;jJ^ and rc=19.4 +2.5 arcseconds. (r,. = 145 + IS/ig^ 
kpc, Table|5}. 

From a spectral extraction region with a radius of 120" and 
using the in-field background subtraction technique, we mea- 
sured the following temperature, metal abundance and redshift; 

= 2.6 + 0.35 keV, Z = 0.38 +0.09, z = 0.57+0.3 (Figure[TT}- 
The overall of the spectral fit was 193 / 214 degrees of free- 
dom. We note the consistency between the redshift measured 
from the X-ray spectrum with the optically determined value 
(z = 0.577; Burke et al., 200 3V The hydrogen column density 
was fixed at the Galactic value (N// - 0.99 x 10^" atoms cm"^), 
when left as a free parameter, its best fit value was very similar 
(Nh = 8.7+^-^ X 10''' atoms cm^^) derived values for the 



bolometric luminosity and absorbed flux inside the viral radius 
are given in Table|5l 

4.4. RXJ0505.3-2849 

The RXJ0505. 3-2849 field is shown in Figurel^d). Four point 
sources removed during analysis are circled and their positions 
listed in Table \K1] Tentative identifications of these sources 
with objects in the UK APM survey ( Maddo x et al., 1990| l sug- 
gest an absolute astrometric accuracy for our observations of 
~2 arcseconds. Figure [T9l shows the corresponding radial sur- 
face brightness distribution and best fit /3 model; yS=0.66;^Q Q4 
and rc=22.8 +2.4 arcseconds. (r^ = 164 + 17/1^^ kpc, Table|5}- 
From a spectral extraction region with a radius of 120" and 
using the in-field background subtraction technique, we mea- 
sured; = 2.5 ± 0.3 keV, Z = 0.17 + 0.08 (Figure[l3. The 
overall Xv of the spectral fit was 279 / 248 degrees of free- 
dom. The fixed Galactic hydrogen value was Nh = 1.5 x 10^" 
atoms cm"^), compared with a free fit parameter, of Nh - 
1.1 + 0.6 X lO^" atoms cm ). The redshift was fixed at its op- 
tically derived value of z = 0.5 1 (j^urke et al ., 2003| l, but when 
left as a free parameter, its best fit value was z = 0.53+0.04. As 
a test, we have also performed a spectral fit without excluding 
the four point sources. We find that, within the errors, the fitted 
temperature was unchanged. 

4.5. RXJ0847.2+3449 (and XMMU 
J084701. 8+34511 7) 

FigurelTja) shows the image of RXJ0847.2+3449. Of the eight 
clusters in our sample, RXJ0847. 2+3449 is noteworthy as be- 
ing the only one that is not a SHARC cluster; it was selected in- 
stead from the 160SD (Vikhl inin et al., 1998t . RXJ0847.2 was 
observed during a period when the PN camera was temporarily 
disabled by a hardware fault, and the requested exposure dura- 
tion was partly compensated by an extension of the MOS obser- 
vation. To improve the contrast of the cluster against the nearby 
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bright QSO (PG 0844+34), only the inner part of the image, 
that covering central CCD of the MOS cameras, is shown in 
Figure Ela). We note that the mirror scattering is sufficiently 
low that no flux from PG 0844+34 impedes our analysis. The 
bright point source ~3 arcminutes west of the RXJ0847.2 clus- 
ter in Figure 0a) was identified with a second, fainter, quasar, 
QSO 0843+349 at z - 1.57. The image also shows evidence 
for a nearby /fl/nfer extended object to the NW, which we iden- 
tify as a cluster candidate with provisional designation XMMU 
J084701. 8+345 117. Vikhlinin (private communication) claims 
to have found a concentration of faint galaxies in this region, 
but to date no galaxy redshifts are available. Neither QSO 
0843+349 nor XMMU J08470 1 .8 +345 1 1 7 fafl inside our spec- 
tral extraction region. The source detection software picked out 
four sources in that region, none of which have counterparts in 
the NED catalogs. These sources were excluded from the spa- 
tial and spectral analysis, they are circled in Fig. and their 
positions are listed in Table lA.fl 

Figure |20l shows the radial surface brightness distribution 
for RXJ0847.2 and the best fit /? model; yS=0.8 1+0.07 and 
Vc=A2 ^4 arcseconds. (r^ - 307 ± 2>Qh^^^ kpc. Table |5}- From 
a spectral extraction region with a radius of 120" and using 
the in-field background subtraction technique, we measured 
Jx = 3.62+°-58 keV, Z = 0.30 ± 0.28 (Figure[T3. The overafl Xv 
of the spectral fit was 152 / 180 degrees of freedom. The fixed 
hydrogen column density was N// - 3.2 x 10^° atoms cm"^, as 
a free parameter, its best fit value was Nh - 2.8 ± 0.14 x 10^" 
atoms cm~^. The redshift was fixed at its optically derived value 
of z = 0.56 but its best fit value was z - 0.54 ± 0.04. As a 
test, we have also performed a spectral fit without excluding 
the four point sources. Doing so changed the fitted temperature 
by 0.1 keV and the measured flux by <10%. 

We have also estimated the temperature of the second clus- 
ter in the field, XMMU J084701. 8+345 117. For this we used 
a spectral extraction radius of 80 arcsec and the in field back- 
ground subtraction technique. We do not have an independent 
estimate of the cluster redshift, and it cannot be constrained by 
the X-ray spectrum, so we have assumed that the cluster lies 
at the same redshift as RXJ0847.2 (z = 0.56). We have also 
fixed the metal abundance to be Z = 0.3 and the hydrogen col- 
umn density to the same Galactic value. We measure a cluster 
temperature of - 1.8^q4 keV and a corresponding (r < r,,, 
0.5-2 keV) flux and bolometric luminosity of 1.3 x 10"'^ 
erg cm"^ s"' and 1 x 10'*-'erg s"' respectively. These results 
are consistent with measured low red-shift - relations. 

4.6. RXJ1 325.5-3826 (and RXJ1 325.0-381 4) 

The data for this cluster were obtained by kind permission of 
Guaranteed Time holder M Watson, from an observation of the 
nearby object IRAS 13224-3809. Unfortunately for our own 
purposes, this bright source had been observed in the MOS 
cameras in a "Window" mode of readout on the central CCDs 
in order to minimize effects of photon pile-up, and conse- 
quently the area of focal plane containing the SHARC cluster 
was not exposed. Therefore, our spatial and spectral analysis 
is restricted to the PN camera data only. Figure0b) shows the 



r 




^^^^^^^^^^^^^^^^^^^^ 




4*1 





0.5 1 2 5 

Energy (keV) 



Fig. 13. Spectral fit and residuals for RXJ0847.2 (MOS data 
only) 



vignetting corrected, background subtracted 0.3^.5 keV PN 
image of the region surrounding IRAS 13224-3809. Clusters 
RXJ1325.5-3826 and RXJ1325.0-3814 (see below) are la- 
beled, as is the streak corresponding to the "out of time" events 
from IRAS 13224-3809. IRAS 13224-3809 is bright enough 
that the "spider" like structure in the PSF is visible around it. 
Given the proximity of the cluster to the IRAS source, the sur- 
face brightness profile for RXJ1325. 5-3826 (Figure l20l i was 
extracted from manually selected clean regions that avoided the 
IRAS 13224-3809 PSF 

For the /5 model fitting, a PSF appropriate for the PN at the 
off-axis angle of the cluster was used (for all the other clusters 
in our sample, we had used the on-axis PSF). The best fit pa- 
rameters in the spatial fit are y6 = 0.64+o q7 and ft. - 17.3 ± 3.0 
arcseconds (r^ = 115 + 2Qh^^ kpc, Table|5j- There is evidence 
for possible excess brightness in the inner bins of the radial pro- 
file. This may be in an indication of a central cooling core, but 
we do not have adequate signal to noise to confirm this spec- 
troscopically. Following the procedure adopted for RXJ1701.3 
(see below) we also performed ayS fit after exclusion of the core 
120kpc and found a value foryS of 0.71+0.05, given a fixed r^ 
of 200/; jQ kpc (We fixed this value arbitrarily to match typical 
cluster values). 

We used a 90" radius extraction region to fit the spectrum 
from this cluster. This is a smaller region than was used for the 
other clusters (120" or 145", Table|5} because of the proxim- 
ity of IRAS 13224-3809. Even though the proton background 
was low enough in this exposure to permit spectral measure- 
ments using the double subtraction technique, we decided to 
use instead only the in-field background subtraction technique. 
This was because we wanted to mimic the off-axis angle de- 
pendence of the point source contamination. Rather than using 
an annulus around the cluster to determine the background, we 
used two source and streak free circular regions at the same 
off-axis angle as the cluster From the background subtracted 
PN spectrum, we measured the following temperature, metal 
abundance and redshift; = 4.15+°-^ keV, Z = 0.31!° ]^, 
z - 0.44 + 0.01 (Figure fT4l . We note the consistency between 
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Fig. 14. Spectral fit and residuals for RXJ1325.5 (PN data 
only) 



the redshift measured from the X-ray spectrum with the op- 
tically determined value (z = 0.445; Burke et al., 2003 1. The 
overall Xv of the spectral fit was 229 / 242 degrees of free- 
dom. The fixed and free hydrogen column density values were 
d(NH - 4.8 X 10^" atoms cm"^ and 4.5+Qg x 10^° atoms cm"^ 
respectively. 

4.6.1. RXJ1 325.0-381 4 

For completeness, we note that another, fainter. Southern 
SHARC cluster lies within the same field of view. The PN im- 
age of this cluster, RXJ1325 .5-3814, can be seen in Figure0b). 
The cluster is also visible in the MOS data (not shown) because 
it lies outside the region blanked by the "Window" mode. The 
cluster is too faint and too far off axis (i.e. where the PSF is 
poorly defined) to permit a detailed spatial analysis, however 
we were able to make a spectral analysis. For this we gener- 
ated off-axis response matrices for each camera using the SAS 
ARFGEN package. From a spectral extraction region with a 
radius of 90", and using the in-field background subtraction 
technique, we measured the following temperature and redshift 
(N// fixed at the same galactic value of 4.8x10^" atoms cm"^ 
and abundance at 0.3) - 3.2 + 0.4 keV and a redshift of 
z - 0.29 + 0.02. We note the consistency between the red- 
shift measured from the X-ray spectrum with the optically de- 
termined value (z - 0.296; Burke et al., 2003 1. 

The cluster is observed close to the edge of the outer CCDs, 
preventing the use of an extraction region large enough to cor- 
respond with the virial radius. Within a radius ~500 /i^^kpc, the 
measured flux and bolometric luminosity are 6.1 (±0.3) 10 ''^ 
erg cm-2 g-' (0.5-2 keV) and = 1.1 + 0.23 lO'^'^gj.g g-i j.g. 
spectively. This combination of Tx and are not consistent 
with the measured - relation for nearby clusters. For ex- 
ample, for Tx = 3.2 keV, one would expect a bolometric lumi- 
nosity of Lx = 2 x lO'^'^erg s"' based on the Markevitch (1998) 
relation, and is probably due mainly to the missing flux outside 
our spectral extraction radius. 



4.7. RXJ1701.3+6414 (and Abell 2246) 

This cluster was observed during an exposure scheduled near 
the end of an orbit, where the spacecraft was approaching the 
particle belts. Not only were the spacecraft operations termi- 
nated prematurely, but most of the exposure was dominated by 
high soft proton background from the edge of the particle belts. 
Fortunately this is one of the brightest clusters in our sample 
and we were still able to produce acceptable quality images 
and spectra. 

Figure ITJc) shows the image of the RXJ 170 1.3 +64 14 
field. Two nearby sources are marked on the Figure; a QSO 
(HS1700-H6416) and a known cluster (A2246). 

The XMM-Newton data indicate an anomalously low r^ 
value in the yS-fit, and this is supported by CHANDRA data 
( IVikhlinin et al. , 2002*) where this cluster was observed with 
a suspected cooling flow central brightness enhancement. We 
therefore excluded the central bins (120 h^^kpc) and fixed the 
rc to the CHANDRA value of 0.5 arcmin (204kpc). Figurel2Ql 
shows the corresponding radial surface brightness distribution 
and resulting best fit j3 model; (3 - 0.64 + 0.05, 0c - 30 arc- 
seconds. This value for /3 is consistent with the CHANDRA 
measurement (0.62 + 0.03). 

From a spectral extraction region with a radius of 120", 
we measured the following temperature, metal abundance and 
redshift (without core excision); Tx = 4.5^[ q keV, Z - 0.24 + 
0.2, z = 0.44+0.02 (Figureinj. The overall Xv of the spectral fit 
was 58 / 53 degrees of freedom. The hydrogen column density 
was fixed at the Galactic value [Nh = 2.6 x 10^*^ atoms cm"^), 
and as a free parameter, its best fit value was very similar (Nh = 
2.6 + 1.5 x 10^° atoms cm"^). 

After excluding the 120kpc cooling flow region we obtain 
Tx = 4.8^[ 3 keV. The derived values for the bolometric lu- 
minosity and absorbed flux inside the viral radius are given in 
Table |5] We were unable to apply the double subtraction tech- 
nique to these data because of the enhanced proton background 
during the observation (see Figure |5jl. The poorer quality of 
data in this observation limits the quality of determination of 
T and L. Agreement with the CHANDRA data is acceptable 
for the temperature, but not the luminosity (5.8+0.5 keV and 
15.7x1 Q*'^eTg/s respectively) . 

4.7.1. Abell 2246 

For completeness we also offer our interpretation of the Abell 
2246 spectrum. From a spectral extraction region with a ra- 
dius of 90", and using the in-field background subtraction tech- 
nique, we measured the following temperature, metal abun- 
dance and redshift by fitting to the PN, MOSl and MOS2 
data simultaneously; = ^-^los ^^^^ ^ = ^-^^ + 0.13, 
z - 0.22 + 0.04. We note the consistency between the red- 
shift measured from the X-ray spectrum with the optically de- 
termined value (z - 0.225; iStruble and Rood, 1991|) . Within a 
radius ~415 /i^g kpc, the flux and bolometric luminosity are 
3.4+0.2 10-1^ erg cm^^ s"' (0.5-2keV) and 2.1 lO^^erg s"' re- 
spectively. Despite the rather restricted radius from which the 
luminosity is determined, the local Lx - relation would pre- 
dict a slightly lower luminosity than this estimate. 
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Spectral fit and residuals for RXJ1701.3; black - PN, 
light grey - MOS 
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Fig. 16. Spectral fit and residuals for RXJ1354.2 black - PN, 
dark & light grey - MOS 



4.8. RXJ1 354.2-0222 

The RXJ1354.2-0222 field is shown in Figure ITJd). A point 
source that was removed during analysis is circled and its po- 
sition listed in Table lA.ll Figure |20l shows the correspond- 
ing radial surface brightness distribution and best fit /3 model 
that gives j6=0.68 ±0.06 and rc=33.6 ^^'^ arcseconds. (rc = 
248!^^/i5^ kpc, Table|5}. 

From a spectral extraction region with a radius of 120", and 
using the in-field background subtraction technique, we mea- 
sured; = 3. 66:^°-^ keV, Z = 0.25 + 0.14, (Figure [16). The 
overall of the spectral fit was 120 / 147 degrees of freedom. 
The fixed Galactic hydrogen column value was N/^ = 3.4 x 10^** 
atoms cm"^, but Nff = 3.2 + 1.2 x 10^" atoms cm"^) when 
free. The redshift best fit value was 0.53 + 0.04. Similarly to 
RXJ1701.3, we were unable to apply the double-background 
subtraction technique to these data because of the enhanced 
proton background (see Figure|5}. 




5 4 5 
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Fig. 17. Lx - relation for an EdS model (^o = 0.5) and 
//()=50 km s ' Mpc ^ The solid line is the best fit for our 
sample after cooling flow excision, the dotted line is the rela- 
tion from Markevitch (1998) for low red-shift clusters. 



5. Discussion 

5.1. Luminosity-Temperature Relations 

Table |5] summarizes our analysis of the eight clus- 
ter observations described above. The mean ele- 
mental abundance (Z = 0.28 + 0.08), core radius 
(r, = 195 + 50/!-^ kpc) and yS ( /? = 0.70 + 0.05), are 
consistent with typical values measured at low redshift (e.g. 
IFukazawa et al., 1998|Jones & Forman, 1999|Mohr, Mathiesen, & Evrard, 19' 
The clusters were derived from well understood surveys 
(SHARC and 160SD) and thus should be representative of 
the cluster population as a whole at these redshifts. We have 
used our observations to determine the Lx - relation at 
0.45 < z < 0.62 and to investigate evolutionary effects. We 
plot luminosity versus temperature, with one sigma errors, for 
the 8 cluster targets in Figure ^\ using the values T'^-^, L'^^ 
of Table 5. The errors are clearly dominated by those of the 
temperature measurements. 

Characterizing the - relation as 



(1) 



and assuming an EdS cosmology, we find a-2.7 +0.4, and 
Le = 15.9^55 X lO'^'^erg s"' using the bisector variant 
of the BCES fitting package ( |Akritas & Bershady, 1996). 
This relation is shown as the solid line on Figure I17l 
The slope of the relation is similar to most previous mea- 
surements, e.g. for an EdS cosmology, a = 2.64 ± 0.27 
( Markevitch, 1998> ; a = 2.33 ± 0.43 (jAllen & Fabian, 1998> ; 
a = 2.88 + 0.15 (|Arnaud& Evrard, (1999)t ; a = 
2.47 + 0.14 (Ikebe et al., 2002); a ~ l82 + 0.32 
( pVfovicki, Soriiig, & Henry, 2002^ . 

To characterize the possible evolution with z, we assume 



(2) 
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Adopting the low redshift - relation from Markevitch 
(1998) {a = 2.64 + 0.27, U = 12.44 + 1.08 x lO^^erg s"' ) 
we performed a minimization on our data points corrected 
by a (1 + z)'* term. The Markevitch (1998) relation was de- 
rived from ROSAT and ASCA observations of 30 clusters at 
0.04 < z < 0.1 (z = 0.05) and is plotted as a dotted Hne 
in Figure On average, the XMM-Newton data points lie 
away from the Markevitch (1998) relation, suggesting an evo- 
lutionary effect whereby clusters of the same temperature were 
more luminous in the past. The best fit value for A is 0.68+0.26 
iXred-^-'^)- We have also determined the BCES best fit for our 
sample when all the points are corrected with this evolution 
to a single redshift (low redshift z - 0.05 of Markevitch) bin, 
and find a = 2.71+0.22, and U = I2.5^_\l x lO'^'^erg s"' . We 
note that we have chosen to compare with Markevitch (1998), 
rather than with the recent study of 82 clusters with ASCA 
temperatures by Ikebe et al. (2002), because the Markevitch 
(1998) method to derive cluster temperatures is closer to our 
own (Section r5.2.3> . thus we try to minimise any systematic ef- 
fects from different data treatment. In Figure^]we re-plot our 
data and the Markevitch (1998) relation after converting them 
to a common cosmology of //q - 70 km s"' Mpc"',QA = 
0.7, Qm = 0.3. In this cosmology we find A=1.52+°-^^. 

For comparison, we include the data points obtained by 
Vikhlinin et al. (2002) from a compilation of Chandra obser- 
vations of high-z clusters (0.39 < z < 1.26) on Figure [Tslfafter 
adjusting to Ho - 70 km s"' Mpc"'). From this Figure, it is 
clear that both the XMM-Newton and Chandra points lie away 
from the Markevitch (1998) relation and that there is no obvi- 
ous systematic off'set between the XMM and Chandra data (see 
Section |5T3}. VikhHnin et al. (2002) also used Markevitch 
( 1 998) as a low redshift benchmark to investigate - evolu- 
tion. From their data, Vikhlinin et al. (2002) found A = 0.6+0.3 
for Qa=0, Qm=1 and a = 1.5 + 0.3 for a Qa = 0.7, Qm = 0.3 
cosmology. 

The value of A has important implications for our under- 
standing of structure formation and cluster evolution. In a self- 
similar model (Kaiser, 1981) cluster X-ray properties are driven 
by gravitational processes, such as shock heating, and a value 
of A~1.5 would be expected in most cosmological models. 
A self-similar EdS model predicts exactly 1.5, and the value 
is only slightly different in a low density concordance model, 
the correction due to the redshift dependence of the virial den- 
sity (Bryan and Norman, 1998/ being small. In a model where 
the cluster X-ray properties are influenced by non-gravitational 
processes, such as energy injection by AGN's or supernovae, 
lower values of A are predicted (e.g. Tozzi & Norman, 2001 ). 
The evolution deduced from our observations assuming an EdS 
model is significantly below the predicted value. Adopting a 
low density concordance model, on the other hand, leads to 
a value consistent with predictions. We conclude that there is 
now evidence from both XMM-Newton and Chandra for an 
evolutionary trend in the Lx - T^x relation. 

Previous studies of Lx - evolution based on either ASCA 
(e.g. [M ushotzky & Scharf, 1997) |AUen & Fabian, 1998| 
INovicki, Sornig, & Henry, 2002, Sadat et al., 1998), 
ROSAT PSPC (e.g . [Fairley et al., 2000) o r Chandra (e.g. 
[Borgani et al., 2001] [Holden et al., 2002T l temperature 
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Fig. 18. U - Tx relation for //o=70 Qa = 0.7,Qm=0.3. The 
dotted horizontal lines are the CHANDRA sample reported 
by IVikhUnin et al., 2002| and the dashed diagonal line the 
Markevitch (1998) Lx - L'x relation corrected to this cosmol- 
ogy. The solid crosses indicate the results from the clusters in 
our study, corrected for this cosmological model, and the diag- 
onal solid line is the best fit for our data. 



measurements found results consistent with no evolution 



( jSadat et al.,199 8 1. When these studies quote values for A, 
these values are generally smaller than those measured by 
us or Vikhlinin et al. (2002), but are usually still at least one 
sigma away from A = 0. Sadat et al. (1998) found a positive 
evolution A = 0.5 + 0.3 for Qm = 1,Qa - and Novicki, 
Sornig, & Henry (2002) A = 1.1 + 1.1 for Qm = 1, Qa = 0, 
and A = 2.1 + 1.0 for Qa=0.7, Qm=0.3 respectively. It is 
clear that we are beginning to probe evidence that Lx - 
evolution, although many more clusters need to be studied, 
and systematic biases (see below) examined in detail, before 
solid conclusions regarding structure formation models can be 
drawn. 



5.2. Data Treatment and Systematics 

We have made every effort to ensure that the Lx - rela- 
tion presented above is robust. We have adopted a uniform ap- 
proach to the spatial and spectral analysis of the eight clus- 
ters in the sample. We have excised regions with clear signs of 
point source contamination before extracting spectra and sur- 
face brightness profiles. We have not used yS-models to calcu- 
late total cluster count rates, but have rather summed up the 
counts within a viral radius. We have tried to ensure that our 
approach to background subtraction does not bias the mea- 
sured cluster parameters (see below). We have used measured 
temperatures to make conversions between count rates and 
fluxes and, where possible, we have used data from all three 
EPIC cameras to derive Lx and values. For RXJ0847.2 and 
RXJ1325.5, data were not available from all three cameras. 
However, we are confident that the derived quantities for these 
two clusters are robust, as we have shown, using the RXJ 1 1 20. 1 
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Fig. 19. Surface brightness distribution of clusters, compared with best y6-model fit (dotted line). Reading clockwise from top 
left: RXJ1120.1,RXJ1334.3,RXJ0337.7 and RXJ0505.3. 



observation, that PN-only or MOS-only spectra are in excellent 
agreement. 

The background subtraction has been complicated by the 
extended nature of the targets and the high proton flare con- 
tamination in some of the observations. We used background 
template files during the spatial analysis. These files covered a 
limited energy range (0.3^.5 keV), to minimize any calibra- 
tion uncertainties. They were filtered using the same rate fil- 
tering criteria as their respective cluster observation and were 
normalized to it using the particle background rates in the CCD 
regions outside the FOV. For various reasons (see Section l33l . 
we did not use the background templates for the spectral anal- 
ysis. Instead, we adopted an in-field background subtraction 
technique. With this method comes the concern that the proton 
background can be over- vignetted. We have checked for this by 
applying a different method, the double subtraction technique, 
to the RXJ1120.1 and RXJ1334.3 observations. We found that 
the results from both techniques agree within the statistical er- 
rors. For the two observations with the worst proton flare con- 
tamination, RXJ1354.2 and RXJ1701.3, we applied a 2% scal- 
ing to the in-field background spectrum to compensate for po- 
tential over- vignetting of the proton background. In Figure |5l 



we demonstrate that this over-vignetting should not present a 
problem for the other clusters in our sample. 

Despite all these quality controls, we still cannot rule out 
the possibility of there being some systematic bias in our re- 
sults. We discuss possible sources of such bias and the potential 
impact on our claimed - evolution below. 

5.2.1 . Impact of Cooling Core Clusters 

Cooling cores have the joint effect of boosting the cluster lu- 
minosity and of reducing the measured global cluster temper- 
ature. These effects have been shown to introduce scatter in 
the Lx - relation and to modify its slope (Markevit ch, 1998] 
[Alien & Fabian, 199 8 1. Various approaches have been taken to 
alleviate the effects of cooling cores on measured Lx - rela- 
tions, e.g. by excluding cooling core clusters from the analysis 
Arnaud & Evrard, (1999), by fitting two temperature models 
( !keb^^r^002r , by masking the central regions from obvi- 
ous cooling core clusters ( |Vikhlinin et al., 2002 1, or by mask- 
ing the central regions from all clusters ( Markevitch, 1998. 1. 
Based on the images (Figures |6l & Q ™d radial profiles 
(Figures[^&|20|) of the eight clusters in our sample, we do not 
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Fig. 20. Surface brightness distribution of clusters, compared with best y6-model fit (dotted Hne). Reading clockwise from top 
left: , RXJ0847.2,RXJ1325.5, RXJ1701.1 and RXJ1354.2. 



expect cooling cores to have much, if any, impact on the mea- 
sured Lx and Jx values; only RXJ1701.3 and RXJ1325.5 show 
evidence of a central surface brightness excess. We have inves- 
tigated this by comparing the results before and after excising 
a circular region * from the cluster cores. 

The results of the re-analysis can be found in Table|5l (rows 
4 and 5). Temperatures changed very little - always within the 
statistical error - after we excised the central r = 50/igQ kpc 
region. In two cases, ( RXJ1120.1 and RXJ1701.3), the lu- 
minosity changes by more than Icr, but still by only 5% and 
9% respectively. Without the excision the best fit - re- 
lation is Le = 17.6fj-3 x lO^^erg s"' and q'=2.66+0.25 for 
an Qa = 0.0, Qm - 1-0, i/o = 50 km s ' Mpc ' cosmol- 
ogy at a typical redshift z ~ 0.55. In summary, we find the 
Lk - Tx relations before and after the core excision are consis- 
tent within the errors. However, we caution that some residual 
cooling flow contamination may still remain after this excision; 
SQh^^^kpc ^ 14", corresponding to an encircled energy fraction 
of a little under 50% of the on-axis mirror PSFs. 

The radius of the excised region was 50 h'^g kpc except in the case 
of RXJ1701.3 and RXJ1325 where ~120kpc was used, see sections 
3.6 and 3.7 



5.2.2. Incompleteness and Flux Errors in the ROSAT 
Catalogs 



Our Lx - evolution result implies that high redshift clus- 
ters in our sample are more luminous than clusters of the same 
temperature at lower redshift. We interpret this as evidence 
for a general evolutionary trend in the Lx - T^x relation, but 
it might also reflect an underlying selection bias in the cluster 
catalogs from which the sample was drawn, e.g. if the SHARC 
and 160SD surveys were biased toward merging systems. It is 
possible that the luminosity and temperature of clusters can 
change dramatically during the process of a major merging 
event (Ricker & Sarazin, 2001,) . Based on the XMM-Newton 
images of the 8 clusters studied here (Figures |6l & 0, such 
a bias does not appear to exist; only RXJ1354 has evidence 
for a possible double component, and excluding this from the 
best fit Lx - ^x relation changes the luminosity normaUsation 
only ~2%. As long as the angular resolution of observations 
allows sufficient discrimination against obvious mergers, then 
presumably an ensemble sample of local and distant clusters 
are similarly affected by merger boosting of luminosity. 
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Moreover, we find no evidence for any incompleteness in 
the SHARC and 160SD surveys, at least in the small area cov- 
ered by these observations. No new clusters were found in the 
fields surrounding the seven SHARC clusters in our sample. In 
the field of the 160SD cluster, RXJ0847.2, we have discovered 
a new X-ray cluster (XMMU J084701.8-H345117). However 
this cluster was too faint to have been included in the origi- 
nal 160SD survey. We estimate the (r < r,,) flux of this object 
to be 1.3 X 10"'"* erg cm"^ s"^ (based on a fitted temperature of 
Tx =1.8 keV and an assumed redshift of z = 0.56), compared 
to the 160SD flux Hmit of 3.7 x lO"'"* erg cm"^ s"' , so this 
cluster is too faint to have been included in the original 160SD 
survey. In summary we do not find evidence for any intrinsic 
bias in the SHARC or 160SD survey selection functions that 
could mimic Lx - evolution. 

Another possible way to mimic the claimed evolution 
would be to systematically overestimate the Lx values. We 
have investigated this by comparing our flux measurements 
to those published in several works ( |Romer et al., 2000| 
IBurke et al., 2003| |\^linin et al., 1998i l. We report here that 
a systematic offset does appear to exist, but in the opposite di- 
rection to that needed to mimic evolution. The XMM-Newton 
fluxes quoted in Table |5l are all lower than previously deter- 
mined values, typically by 20% but as much as by 40% in the 
case of RXJ0505.3. We note that this discrepancy persists re- 
gardless of the ROSAT catalog from which the cluster was se- 
lected. The average discrepancy between the Southern SHARC 
(4 clusters). Bright SHARC (3 clusters) and 160SD (3 clusters) 
fluxes and the XMM-Newton fluxes is 22%, 21% and 16% re- 
spectively. 

Improvements in point source rejection, count rate to flux 
conversions and the signal to noise have meant that the typ- 
ical statistical error on the XMM-Newton determined fluxes 
in Table m is ^ 3%, compared to ^ 10% for the ROSAT de- 
terminations. Therefore, even accounting for a possible offset 
in the absolute cross calibration of the two observatories, it is 
clear that a systematic flux offset does exist. There are likely 
to be several reasons for this offset. Of these, the use of an 
incorrect spatial model, to convert between an aperture flux 
and a total flux, is probably the most significant. As shown in 
IRomer et al., 2000, the use of a generic (j3 - 0.67, rc - 250/!g(J 
kpc), rather than best fit, y6-model resulted in a typical enhance- 
ment of 10% in the total cluster flux. In our XMM-Newton 
analysis, we do not use the j6-model fits to determine the total 
flux. Instead we directly sum up the counts within a virial ra- 
dius. For some clusters this difference in approach can explain 
the entire ROSAT to XMM-Newton discrepancy. For example, 
in the case of RXJ1334.3, our estimate of r,, is almost identi- 
cal to the rgo radius within which 80% of the total flux from 
a generic j6-model would fall. Romer et al., 2000. divided the 
ROSAT counts within rgo by 0.8 to estimate the total count rate. 
For RXJ 1334.3 this division would artificially boost the flux by 
25% which more than accounts for the 22% the discrepancy be- 
tween die XMM-Newton and Bright SHARC values. 

For the four Southern SHARC clusters in our sample, we 
have recalculated the total ROSAT PSPC fluxes using the /? 
and rc values given in Table |5l (in |Burke et al., 2003| generic 
values were used; (3 - 0.67, rc - 250h^Q kpc). By doing so. 



we improve the consistency between the ROSAT and XMM- 
Newton fluxes to within 5% for RXJ 1354.2 and RXJ 1325 .5 
and within 15% for RXJ0337.7. For RXJ0505.3 there remains 
a large, 35%, discrepancy which requires further investigation. 

For the three Bright SHARC clusters in our sample, we 
have gone back to the original Romer et al., 2000| data and mea- 
sured count rates inside r,,. Doing so reduces the ROSAT de- 
termined flux by ^ 30% in the case of both RXJ 1120.1 and 
RXJ1334.3, i.e. so that they are ^ 10% lower than the val- 
ues quoted in Table [S] In the case of RXJ 170 1.3, the measured 
ROSAT flux actually increased slightly when measured inside 
ry rather than rgo. On closer examination of the ROSAT data, 
it was apparent that the nearby QSO and cluster (see FigureQ 
were contaminating both the source and background apertures 
in the PSPC image. When the source and background were ac- 
cumulated only from the (source free) region to the West of the 
cluster, the ROSAT flux within r,, dropped to within ^ 10% of 
the XMM-Newton determined value. These examples demon- 
strate the importance of using XMM-Newton to re-calculate 
fluxes and luminosities for high redshift clusters detected at 
low signal to noise by ROSAT. We stress that, despite these flux 
uncertainties, the SHARC and 160SD catalogs are still fair rep- 
resentations of the underlying cluster population and can still 
be used to probe the - relation. It is also important to note 
that the these flux uncertainties do not apply to the ROSAT ob- 
servations used by iMarkevitch, 1998| to determine Lx values for 
the 30 low redshift clusters in his Lx - Ty_ analysis. Those ob- 
servations have exquisite signal to noise, so that the Lx values 
derived from them will be limited only by the absolute calibra- 
tion of the instrument. 

The fact that we find a Lx - Lx relation consistent with the 
low redshift value argues against a bias in our selection towards 
the more luminous tail of the cluster population; for either 

- all the clusters must be biased in the same way, or 

- they are biased in just such a way as to offset any real evo- 
lution in the value of the slope. 

Nevertheless the reader might be cautioned that should 
there remain a common bias for all cluster samples, despite 
these cross-checks, it is possible the clusters selected represent 
the brightest portion of the intrinsic high-z sample, and thus re- 
main on the high luminosity end of the true population, hence 
mimicking the luminosity evolution we see. Further valida- 
tion of the ROSAT data with CHANDRA and XMM-Newton 
should help to close this issue. 

5.2.3. Comparison with IVIarkevitch (1998) and 
Vikhlinin et al. (2002) 

We come to very much the same conclusions regarding Lx - Lx 
evolution as did |Vikhlinin et al., 2002| Given that these con- 
clusions provide evidence for significant evolution, contrary to 
previous analyses, it is important to investigate whether this 
concordance is genuine or coincidental. Both studies rely on 
'Markevitc h, 1998| for the low redshift comparison, so let us first 
consider how our XMM-Newton clusters might appear either 
too luminous, or too cool, compared to the [Markevitch, 1998| 
sample: 
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1. Markevitch 1998 measured cluster fluxes by summing the 
count rate inside circular apertures. However, rather than 
using a cluster specific virial radius, he used a fixed metric 
aperture of r - l/ig^J Mpc and ignored any flux that might 
lie outside this radius. However, this would be a very small 
efl'ect; Vikhlinin et al. (2002) estimate that the no more than 
4% of the total flux would be missed by adopting a fixed 
r = l/igy Mpc aperture. 

2. When Markevitch 1998 excised the central r = SOh^^ kpc 
region from all the clusters in his low redshift sample be- 
fore measuring the flux, he applied a multiplication factor 
of 1.06 to account for flux falling inside the excised re- 
gion. We carried out an excision and flux lost correction on 
our eight clusters with a cluster-specific surface brightness 
model. (Section r5.2.1l and TableO, we found that for six of 
the eight clusters, the measured values dropped slightly. 
The difference in our techniques could systematically affect 
cluster brightness, compared with the low redshift counter- 
parts, but only at the few percent level. 

3. Uncertainties in the cross calibration between XMM- 
Newton and the instruments used for the |Markevitch, 1998| 
analysis (ROSAT and ASCA) may mimic evolution in the 
Lx - Tx relation. At this time it is not possible to rule out 
that possibility; to date, a full comparison of XMM-Newton 
and ASCA determined cluster temperatures has not been 
carried out. 

4. Ikebe et al., 2002 have used a different approach to 
[Markevitch, 1998| to analyse ASCA observations of low 
redshift clusters. They measure temperatures that are on 
average lower than those measured by [Markevitch, 1998] 
with the trend becoming more pronounced as increases. 
However, this should not impact our conclusions regarding 
Lx - evolution, given that we used a similar technique to 
[Markevitch, 1998| to measure I v . 

5. Some estimate of the effect of external systematic error 
could be taken by considering the value of A in fitting to a 
different low-z - sample. We have done so using the 
[Arnaud & Evrard, (1999)| relation and find A = 0.95 ± 0.2 
(again in the Einstein de Sitter cosmology). Therefore, de- 
spite a possible contamination in that sample from cooling 
flows we still see evidence for evolution 

We note that the first two concerns listed above do not 
apply to the Vikhlinin et al., (2002) study, as they deliber- 
ately adopted the same methodology as [Markevitch, 1998| and 
the last three concerns are not likely to be important either; 
Vikhlinin et al., (2002) have carried out a cross comparison of 
Chandra and ASCA values for ^ 20 low redshift clusters and 
find them to be in agreement at better than the 5% level, with no 
systematic offset. Despite this, we should still consider the pos- 
sibility that both studies (ours and [ Vikhlinin et al., 2002[ l might 
be pointing to the wrong value of A. At present we cannot quan- 
tify how factors such as uncertainties in the cross-calibration of 
Chandra and XMM-Newton, especially at low energies, might 
impact our conclusions. We note in particular that Vikhlinin 
et al., (2002) applied an empirical factor of 0.93 to improve 
the cross-calibration of the ACIS front- and back-illuminated 
CCDs below 1.8keV. Other differences between the observa- 



tories may also be important. For example, it will be easier to 
account for point source and cooling core contamination using 
Chandra data, because of the improved spatial resolution (see 
for example the case of RXJ1701.3). Alternatively, it should be 
easier to correctly account for the various particle (cosmic ray 
and proton) backgrounds using XMM-Newton data, because 
its CCD's cover a larger area both inside and outside the FOV 
In summary, although we cannot rule out the possibility that a 
combination of factors have lead us to measure an artificially 
large value for the A parameter, we believe that our XMM mea- 
surements provide clear evidence for evolution in the - 
relation. 
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Table 5. Summary of cluster parameters for EdS model Ho=50, qo=0.5. Spectral fitting errors, L,y6 and r^. are Icr on one 
parameter. Fluxes are the measured, absorbed fluxes in ROSAT band. Fractr is the fractional correction made from the spectral 
extraction region, to the total counts within the r,, radius after accounting for point source excision, inter-chip gaps loss etc.. 
The 3rd & 4th rows summarise the data for the case when the core region has been excised from the cluster core. (Except for 
RXJl 70 1.3, where following the CHANDRA data, we excise a cooling flow enhancement to 120kpc). 
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Appendix A: Excised Point Sources 



XMM-Newton ID 
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8-252318 
7 - 252326 
9 - 252208 
9 - 252206 
2 - 252214 



03:37:37.8 
03:37:43.7 
03:37:42.9 
03:37:45.9 
03:37:47.2 



■25:23:18.0 

■ 25:23:26.7 
■25:22:8.4 

■ 25:22:6.2 
■25:22:14.9 



Field RXJ0505.3 
XMMU J050512, 
XMMU J050517, 
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XMMU J050510, 
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5 -285034 
3 ~ 285023 
3 - 285006 
9 - 284951 



05:05:12.5 
05:05:17.3 
05:05:22.3 
05:05:10.9 



-28:50:34.6 
-28:50:23.7 
- 28:50:06.3 
-28:49:51.1 



Field RXJ0847.2 
XMMU J084711, 
XMMU J084714, 
XMMU J084707, 
XMMU J084709, 



+ 3449 

4 + 344717 

5 + 344654 
5 + 344947 
5 + 344917 



08:47:11.4 34:47:17.1 

08:47:14.5 34:46:54.4 

08:47:7.5 34:49:46.7 

08:47:9.5 34:49:17.1 



Field RXJ 1120.1 
XMMU Jl 11959, 
XMMU Jl 12001. 
XMMU Jl 12004. 
XMMU Jl 12008, 
XMMU Jl 12009, 
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Field RXJ 1334.3 
XMMU J 133426, 
XMMU J133410, 
XMMU J133415, 
XMMU J133416, 
XMMU J133430, 
XMMU J133428, 



+ 5030 
3 + 503247 

8 + 503118 
6 + 503030 
3 +503115 
2 + 503238 

9 + 503141 



13:34:26.3 
13:34:10.8 
13:34:15.6 
13:34:16.3 
13:34:30.2 
13:34:28.9 



50:32:46.8 
50:31:17.6 
50:30:29.7 
50:31:15.4 
50:32:38.0 
50:31:41.5 



Field RXJ1354.2 
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Table A.l. The identifications (XMM-Newton informal ID, nominal RA & Dec ) of point sources that were excised from the 
spectral and imaging analysis 



